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MIXED EHRHART POLYNOMIALS 


CHRISTIAN HAASE, MARTINA JUHNKE-KUBITZKE, RAMAN SANYAL, 
AND THORSTEN THEOBALD 


Abstract. For lattice polytopes Pi,... ,Pk Q Bihan (2016) introduced the discrete 
mixed volume DMV(Pi,..., Pk) in analogy to the classical mixed volume. In this note we 
study the associated mixed Ehrhart polynomial MEpj_..,_p,,(n) = DMV(nPi,..., nPfc)- 
We provide a characterization of all mixed Ehrhart coefficients in terms of the classical 
multivariate Ehrhart polynomial. Bihan (2016) showed that the discrete mixed volume 
is always non-negative. Our investigations yield simpler proofs for certain special cases. 

We also introduce and study the associated mixed h*-vector. We show that for large 
enough dilates rPi,... ,rPk the corresponding mixed h*-polynomial has only real roots 
and as a consequence the mixed h*-vector becomes non-negative. 


1. Introduction 

Given a lattice polytope P C the number of lattice points |P H Z‘^| is the discrete 
volume of P. It is well-known ([11]; see also [2, 3]) that the lattice point enumerator 
Ep{n) = lnPnZ‘^1 agrees with a polynomial, the so-called Ehrhart polynomial Ep{n) G 
Q[n] of P, for all non-negative integers n. Recently, Bihan [ 6 ] introduced the notion of 
the discrete mixed volume of k lattice polytopes Pi,..., P^ C 

(1) DMV(Pi,...,Pfc) := 5^(-l)'=-l^l|PjnZ''|, 

JC[A;] 

where Pj := Yhj&jPj for 0 7 ^ J C [k] and P 0 = {0}. 

In the present paper, we study the behavior of the discrete mixed volume under simul¬ 
taneous dilation of the polytopes Pj. This furnishes the definition of a mixed Ehrhart 
polynomial 

(2) MEp,,..„p,(n) := DMV(nPi,..., nP^) = (-l)'=-l'"Pp,(n) G Q[n]. 

JC[fc] 

Khovanskii [16] relates the evaluation MEp^ . ,_pj,(—1) to the arithmetic genus of a compac- 
tified complete intersection with Newton polytopes Pi,..., Py- Danilov and Khovanskii [9] 
investigate the Hodge-Deligne polynomial e{Z-,u,v) G 7j[u,v] of a complex algebraic va¬ 
riety Z which in case of a smooth projective variety agrees with the Hodge polynomial 
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author verified that DMV(Pi,..., Pk) = 1, 0) where Z C (C*)'^ is the uncom- 

pactified complete intersection [1], Using the mixed Ehrhart polynomial, this yields the 
reciprocity type result 

MEp,,.,„p,(-l) = {-iy-’^eiZ-l,0) while MEp,,..„p,(l) = (-l)"-'^e(Z; 1, 0), 


relating Z and its compactification Z. 

To the best of our knowledge, the origin of the mixed Ehrhart polynomial (2) goes back to 
Steffens and Theobald [22], see also [21], In their work, a slight variant of (2), yet under 
the same name, has been employed as a very specific means to study higher-dimensional, 
mixed versions of Pick’s formula in connection with the combinatorics of intersections of 
tropical hypersurfaces. The definition of the mixed Ehrhart polynomial therein differs 
from our definition by the exclusion of the empty set from the sum. 

Here, we study mixed Ehrhart polynomials and their coefficients with respect to various 
bases of the vector space of polynomials of degree at most d. First, we show that in 
the usual monomial basis, the coefficients of mixed Ehrhart polynomials can be read off 
directly from the multivariate Ehrhart polynomial Pp(ni,...,n^), where ni,... ,nk € Z>o 
are non-negative integers (Theorem 2.4). This gives a meaning to the coefficients of 
MEpj_,.._Pj,(n). In particular this provides a simple proof that the coefficient of n* vanishes 
ioY i < k and also allows to give streamlined proofs for known characterizations of the two 
leading coefficients (Corollaries 2.5 and 2.6). We then deal with two prominent subclasses. 
For the case that all polytopes Pi,... ,Pk are equal, the mixed Ehrhart polynomial can 
be expressed in terms of the h*-vector of P (see Proposition 2.7). And for the case that 
Pi,... ,Pk all contain the origin and satisfy dim(Pi Pk) = dim Pi dimP^, 

we can provide a combinatorial interpretation of DMV(Pi,..., Pk) (see Proposition 2.9). 
As a consequence of the main result from [6], it follows that MEp^ . Pj,(?7,) > 0 for all n > 0 
(see Theorem 3.1). This is accomplished by the use of irrational mixed decompositions 
and skilled estimates. Using our understanding of the coefficients of MEp^., Pj,(n), we 
give direct proofs for this fact in the cases k G {2, d — l,d} and Pi = ■ ■ ■ = P^ = P in 
Section 3. See also [14] for further developments in the context of valuations. 

Expressing the (usual or mixed) Ehrhart polynomial in the basis : 0 < i < d} 

gives rise to the definition of the (usual or mixed) h*-vector and h*-polynomial (see Sec¬ 
tion 4). By a famous result of Stanley [18], the usual h*-vector is non-negative (compo¬ 
nentwise). We illustrate that for the mixed h*-vector this is not true in general (see Exam¬ 
ple 4.3). Yet, we show that this has to hold asymptotically for dilates rPi,rP2 ,... ,rPk 
(Corollary 4.6). This follows from the stronger result that for r ^ 0 the mixed h*- 
polynomial is real-rooted with roots converging to the roots of the d^^ Eulerian polyno¬ 
mial (Theorem 4.5). This can be seen as the mixed analogue of Theorem 5.1 in [10] (see 
also [4, 8]). As a byproduct, we obtain that asymptotically the mixed h*-vector becomes 
log-concave, unimodal and, as mentioned, in particular positive, except for its 0**^ entry, 
which always equals 0 (see Corollary 4.6). 

Our paper is structured as follows. In Section 2, we prove various structural properties of 
the mixed Ehrhart polynomial. In Section 3, we review Bihan’s non-negativity result of the 
discrete mixed volume, particularly from the viewpoint of the mixed Ehrhart polynomial, 
and provide alternative proofs for some special cases. Finally, in Section 4 we study the 
h*-vector and the h*-polynomial of the mixed Ehrhart polynomial, and in particular show 
real-rootedness of the mixed h*-polynomial and positivity of the h*-vector for large dilates 
of lattice polytopes. 
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2. Structure of the mixed Ehrhart polynomial 

In this section, we collect basic properties of the mixed Ehrhart polynomial. For some 
known results we provide new or simplihed proofs. For the whole section, we £x a collec¬ 
tion P = (Pi,..., Pfc) of k lattice polytopes in and we assume that Pi P^ is of 

full dimension d. The mixed Ehrhart polynomial, as introduced in (2), is by dehnition a 
univariate polynomial of degree < d, which can be written as 

MEp(n) = med(P)n'^ + merf_i(P)n‘^“^ -I- ■ ■ ■-I- meo(P). 

Example 2.1. For the case of k copies of the d-dimensional unit cube, Pi = ■■■ = Pk = 
[0, we have Ep^{n) = {n + and thus 

MEp(n) = ^(-lUTOon+l)- = E (^) (jn)' 

j=o VJ/ i=o vv 

by the binomial theorem. Hence, 

MEp(.) = 

which gives mei(P) = (i)'^^/(0) with f{x) = P, where A^f{0) denotes the k^^ difference 
of the function f at 0 (see, e.g., [19, Sect. 1.4]d This can be expressed as mei(P) = 
{^^k\S{i,k) in the Stirling numbers S{i,k) of the second kind (see [19, Prop. 1.4.2]d 

Remark 2.2. If P = (Pi,..., Pk) is a collection of lattice polytopes in and dimP* = 0 
for some i, then MEp(?7,) = 0. To see this, assume that Pk is 0-dimensional. For any 
J C [/c —1], the polytopes Pj and Pju{A:} dxe translates of each other and the corresponding 
terms in (2) occur with different signs. 

As our hrst result, we give a description of the coefficients of the mixed Ehrhart polynomial 
in Theorem 2.4. For this we recall a result independently due to Bernstein and McMullen; 
see [13, Theorem 19.4]. 

Theorem 2.3 (Bernstein-McMullen). For lattice polytopes P = (Pi,... ,Pk) in Mf, the 
function 

Ep{ni,... ,nk) := |(niPiH- \-nkPk)FZ‘^\ 

agrees with a multivariate polynomial for all ni,... ,nk € Zi>o. The degree of Pp in Ui is 
dim Pi . 


In particular, Pp(ni,... ,nk) has total degree d = dim(Pi Pk). We write this 

polynomial as 

(3) Ep{ni,. .. ,nk) = ^ e„(P)n“ 

where rP = ■ ■ ■ nff for a = (oi,..., ak) G Z>g. 

Theorem 2.4. For P = (Pi, ..., Pk) a collection of lattice polytopes in 

mei(P) = y^ea(P), 

a 

where the sum runs over all a G such that |q!| := Oi -f ■ ■ ■ -f a*, = i. In particular, 
mei(P) = 0 for all 0 < i < k. 
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Proof. For a subset J C [k] we write Ij G {0,1}*^ for its characteristic vector. Observe 
that 

DMV(Fi,...,Pfc) = 

JC[k] 

Setting 

k 

^(«) = where = 

JC[k] i=l 

presentation (3) then implies 

DMV(Fi,...,Pfc) = e„(P)Al(«). 

Using that 0*^ = 1, it is easy to verify that A{a) = 1 if Oi > 0 for all i = 1,... ,k and 
A{a) = 0 otherwise. Hence, 

d 

MEp,_p^{n) = DMV{nP,,...,nPk) = eo.{nP)A{a) = Y.iYl ^o.iP))n\ 

*=0 

\a\=i 

where the last step used that is homogeneous of degree (ai,..., ad). □ 

The proof also recovers Lemma 4.10 from [22] with different techniques. In order to give 
an interpretation to some of the coefficients of MEp(n), let us write 

d 

Fp(ni,... ,nfc) = ^Fp(ni,... ,nfc) 
i=0 

where Fp(?7,i,..., Uk) is homogeneous in rii,..., rifc of degree i. By the same reasoning as 
in the case of the usual Ehrhart polynomial (see, e.g., [3, Sect. 3.6]), we note 

Ep[ni..... rifc) ^ vold(niPi H-1- n^Pb) 

\a\=d 

In the above expression, the (normalized) coefficient MV is called the mixed volume; 
see, e.g., [13, Chapter 6]. Moreover, the notation MVrf(Fi[ai],..., Ffc[afc]) means that 
the polytope Pi is taken times, and (4) comes from a well-known property of mixed 
volumes (see [17, Section 5.1]). Based on (4) together with Theorem 2.4, we obtain 
streamlined proofs for the following known characterizations of the two highest mixed 
Ehrhart coefficients (see [22] for the case k = d — 1 and [21, Lemma 3.7 and 3.8]). 

Corollary 2.5 (Steffens, Theobald). Let P = (Fi,..., F^) he lattice polytopes such that 
d = dim(Fi F^). Then 

mej(P) = 

Q; 

where the sum runs over all a G with |a| = d. In particular, MEp(n) is a polynomial 
of degree exactly d. 

Proof. By Theorem 2.4, merf(P) is the sum of the coefficients of Fp(?7,i,... ,nd) of total 
degree d. These are coefficients of Fp(ni,... ,nk) = vol(niFi n^Pk) that are 

given by the mixed volumes. Eor the second claim, it is sufficient to note that since the 
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sum Pi + ■ ■ ■ + Pfc is full-dimensional, Pp is not identically zero and all coefficients are 
non-negative. □ 


For k = d, the previous corollary recovers a result of Bernstein [5] 
(5) MEp,_p^{n) = d\ MVrf(Pi,..., P,) n" 


A similar reasoning also allows us to give an expression for the second highest coefficient 
med-i(P). For an integral linear functional a G the rational subspace C is 

equipped with a lattice and a volume form vol^ so that a fundamental parallelpiped 

has unit volume. We denote the face of a polytope P C where a is maximized by P“ 
and consider it (after translation) as a polytope inside a~^. Similarly, we write MV^ 
(DMVa) for the (discrete) mixed volume of a collection of polytopes in a-*- computed 
using vola. 


Corollary 2.6 (Steffens, Theobald). Let P = (Pi,..., P^) he a collection of d-dimensional 
lattice polytopes in Then 


med-i(P) = 


\a\=d—l 


d-i 


Ol, 


Vv,(Pi“[ai],. 


•,^fcK]), 


where a ranges over the primitive facet normals of Pi + --- + Pk- 


Proof. For a full-dimensional lattice polytope P C the second highest Ehrhart coeffi¬ 
cient can be expressed as 

a 

where a ranges over all primitive vectors of facets of P (see, e.g., [3, Thm. 5.6]). This, of 
course, is a finite sum as P has only finitely many facets. It follows that 

JC[k] a 

a JC[fc] 

Since ed{P) = vo\d{P) for a d-dimensional lattice polytope P C M'^, (2) implies that the 
inner sum in the previous expression equals med_i(Pi“, • • • > ^k)- The result now follows 
from Corollary 2.5 applied to P“ = (Pi“,... ,Pfc)- Observe that if a is a facet normal 
for Pj, then a is a facet normal for Pj for all I J. Hence, the above sum is over all 
primitive facet normals of Pi P*,. □ 


Next, we will show that for the special case Pi = P 2 = ■ ■ ■ = P^ = P, the mixed Ehrhart 
polynomial can be expressed in terms of the h*-vector of P; see Section 4. For now, let us 
mention that Stanley [18] showed that there exist non-negative integers h*{P) such that 

( 6 ) E,(„) = 

In terms of this presentation we get: 

Proposition 2.7. Let P be a d-dimensional lattice polytope and let P = (P,..., P) 
he a collection of k copies of P. Then 
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MEp(n) = 

j=0 \i=0 


k-ifk\ fin + d-j 


In particular, 


DMV(P) = 

j=0 


d 




h’AP). 


Proof. From the definition of the mixed Ehrhart polynomial, we infer 


MEp(n) = ^(-1)"-* . E,p(n). 


i=0 


Using that Eip{n) = Ep{in) and the expression of Ep{n) as in (6), this yields 

d k 


MEp(n) = 5^/!*(/>) 5^(-l) 

j=0 


k_ifk\fin + d- j 


i=0 


d 


This shows the first claim. For the second claim it suffices to observe that DMV(P) = 
MEp(l) and to check that, in this case, the inner sum in the above expression equals 
by the binomial identity in [12, (5.24)]. □ 


This is reminiscent to the relation between /-vectors and h-vectors in the enumerative 
theory of simplicial polytopes. Investigating this analogy further yields a theory of discrete 
mixed valuations; see [14]. 

Example 2.8. Let P = [0,1]^ C be the 3-dimensional unit cube. Then the usual 
h*-vector of P equals (1,4,1,0). By Proposition 2.1, the discrete mixed volume of the 
collection (P, P) is given by 


DMV(P, P) 


3-0 

3-2 


■ 1 + 


3-1 

3-2 


■4 + 


3-2 

3-2 


■ 1 


12 . 


The mixed Ehrhart polynomial of (P, P) equals MEp_p(f) = 6t^ -\- Qt^, which is consistent 
with Example 2.1. 


Proposition 2.9. Let P = (Pi,..., P^) be a collection of lattice polytopes, each containing 
0, such that dim(Pi -!-■■■+ Pk) = dim Pi dimP^. Then DMV(Pi,..., Pk) counts 

the number of lattice points z G (Pi -!-■■■+ Pk) D that are not contained in a subsum 
PjforJ^[k]. 


Proof. By the hypothesis, P = Pi-|- - ■ --l-Pfc is affinely isomorphic to the Cartesian product 
Pi X ■ ■ ■ X Pk. For any J C [fc], this implies that Pj is the intersection of P with the linear 
span of Pj. And for I,J^ [k], we see that the intersection P/flPj is exactly Pmj- Hence 

DMV(Pi,...,Pfc) = IPnz'^l - ^(-i)^-i^i|Pjnz‘^| = |(p\ IJ Pj) nz'^l. □ 

J¥=lk] 

For the case that all polytopes are segments, the previous proposition can be used to 
recover (5). 
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3. Non-NEGATIVITY OF THE DISCRETE MIXED VOLUME AND THE MIXED EHRHART 

POLYNOMIAL 

Bihan [6] showed the following fundamental non-negativity result. 

Theorem 3.1 (Bihan). Let Pi ,C be lattice polytopes. Then 

DMV(Fi,...,Pfc) > 0. 

Thus the mixed Ehrhart polynomial evaluates to non-negative integers for all positive 
integers n. 

For a proof of Theorem 3.1, Bihan develops a theory of irrational mixed decompositions 
which yields a technical induction on dimension. The discrete mixed volume is a particular 
combinatorial mixed valuation in the sense of [14]. It is shown in [14] that the discrete 
mixed volume is monotone with respect to inclusion, which implies Theorem 3.1. The 
proofs in [14] are less technical but set in the context of the polytope algebra. In this 
section we give simple and geometrically sound proofs for the special cases k G {2, d — 1, d} 
and Pi = ■■■ = Pk = P ■ For the case k = d this is a consequence of (5). For Pi = 
■ ■ ■ = Pk = P this follows from Proposition 2.7 together with Stanley’s result on the 
non-negativity of the h*-vector; see [18, 15]. 

Direct proof of Theorem 3.1 for k = 2. Let Pi and P 2 be lattice polytopes in M'’*. Since 
the mixed Fhrhart polynomial is invariant under translation of the polytopes, we may 
assume that 0 is a common vertex of Pi and P 2 , Pi fl P 2 = {0} and that there is a 
hyperplane H weakly separating Pi from P 2 . Hence Pi U P 2 C Pi -|- P 2 . It follows that 

Ep^+p^in) > Ep^{n) + Ep^{n) - 1 , 

for all n G Z>o, i.e., MFp^^P 2 (n) > 0 for all n G Zi>o. □ 

For A; = d — 1, a direct proof idea has already been used in the framework of tropical 
geometry and d-dimensional Pick-type formulas in [22]. The following is a proof along 
similar lines. 

Direct proof of Theorem 3.1 for k = d — 1. By Theorem 2.4, the mixed Ehrhart polyno¬ 
mial of a collection of lattice polytopes P = (Pi,..., Pd-i) P is of the form 

MEp(n) = merf(P)n'^-|-merf_i(P)n'’*“^. 

It suffices to show that merf(P) and merf_i(P) are non-negative. Since the mixed volume 
is multilinear and symmetric in its entries. Corollary 2.5 yields 

me,(P) = fMVrf ^Pi,P 2 ,...,Prf_i,^P.j . 

For the remaining coefficient. Corollary 2.6 together with the fact that k = d — 1 yields 


where the sum is over all primitive facet normals of Pi P^-i. 


□ 
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4. Mixed /i*-polynomials 


For a (i-dimensional lattice polytope P, let 

Elirp(z) = Ep{n)z"' 

n>0 


be the Ehrhart series of P. As Ep{n) is a polynomial in n of degree d, there exist 
integers hl{P ),..., h^(P) such that 


Ehrp(2;) 


Y.Uh](P)^l 


see, e.g., [2, 13]. The polynomial in the numerator is called the h.*-polynomial of P 
(also known as (5-polynomial [20]) and denoted by h*p{z). Similarly, the coefficient vector 
h*{P) := (hQ(P),..., h*^{P)) is called the h.*-vector of P. Stated differently, the Ehrhart 
polynomial of a cf-dimensional lattice polytope P can be expressed as in (6). By Stanley’s 
non-negativity theorem [18], the coefficients hg(P), h^(P),..., h*^{P) are non-negative in¬ 
tegers. Moreover, it is known that 


(7) ho(P) = i, h];(P) = iPnz'^l - (d + 1), h;)(P) = |int(P) 


where int(P) denotes the interior of the polytope P. 

The dehnition of mixed Ehrhart polynomials prompts the notion of a mixed h*-vector 
h*(P) = (/iq(P), ..., h^(P)) of a collection P = (Pi,..., P^) of lattice polytopes, which is 
given by 

(8) MEp(n) = ft;(P) (" ^ '*) + • • • + ftj(P) ([[). 


where d = dim(Pi -|- ■ ■ ■ + Pk). In this section, we study properties of h*(P) for a collection 
of full-dimensional lattice polytopes, i.e., dimP, = d for all i. Though Stanley’s h*-non- 
negativity does not extend to the mixed h*-vector as we will see below, we show that 
for large enough dilates rP := (rPi,..., rP^), r 3> 0, the corresponding mixed h*- 
polynomial 

= A;(rP) + ... + A^(rP)z‘' 

has only real roots and hence the mixed h*-vector is log-concave. This is in line with 
results of Diaconis and Fulman [10] for the usual h*-vector. 


Remark 4.1. Note that the contribution of the index set J = 0 in the discrete mixed 
volume (1) is (—1)^. For the mixed Ehrhart series X]n>o this induces 

for the index set J = 0 the contribution 

(9) ^(-1)L" = = (-F 

n>0 


(1 - ^ ’ (1 - 


Since we assume that all polytopes Pi have the same dimension, linearity and Remark 4.1 
allow to write the mixed h*-vector as 

(10) ft-(P) = 5; + (-i)'-'+‘ 

0^JC[fc] 

for 0 < i < (i. The next lemma collects some elementary properties of mixed h*-vectors. 

Lemma 4.2. Let P = (Pi, P 2 ,..., P^) C 6e a collection of d-dimensional lattice 
polytopes. Then: 

(i) /! 5 (P) = 0 , 
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(ii) If k = 2, then h\[Pi^P2) = DMV(Pi,P 2 )- In particular, h\[Pi,P2) > 0. 

Proof, (i) Since, for any lattice polytope P C we have hQ{P) = 1 (see (7)), it follows 
from (10) that 

a;(p)= E (-i)‘-m + (-i)‘=o, 

07^JC[fc] 

(ii) We know from (7) that for a Pdimensional lattice polytope P C 

h*(P) = |P n Z^\ - (d + 1 ) = Ep{l) -{d+l). 

This fact combined with (10) yields 


h*(Pi,P2) = (Pp,+P,(l) - (d+l)) - (Pp,(l) - (d+ 1) +Pp,(l) -{d+l))-d 
= Epi+P2{I) ~ -^Pi(l) ~ Ep.^{l) + 1 
= DMV(Pi,P2), 

where the last equality follows from (2). For the second claim, it suffices to note that 
DMV(Pi,P 2 ) is non-negative (see Section 3). □ 

A natural question when dealing with integer vectors is whether all entries are non¬ 
negative. The previous lemma provides some positive results in this direction. However, 
as opposed to the h*-vector of a single lattice polytope, the next example shows that in 
general it is not true that the coefficients of the mixed h*-vector are non-negative. 

Example 4.3. Consider the collection P = (A^,..., A^) consisting of k copies of the 
standard d-simplex. For{d+l—i)k < d+l none of the dilates {d+l—i)Pj = (d+l—i)|J|Arf 
has interior lattice points. Hence, h*{Pj) = 0 so that h*(P) = (—1)^+*(^) can be negative 
(see (10)). A specific case is h*(A 3 ,A 3 ) = (0,3,4,—1). 

Though we have just seen that there do exist collections of polytopes with negative 
mixed h*-vector entries, observe that for m > 2 all entries of h*{m/S.^,mA^) = (0,m^ + 
2m,4m^,m^ — 2m?) are indeed non-negative, and the leading coefficients in m are the 
Eulerian numbers (0,1,4,1) (see below). 

We therefore propose to study the following question. 

Question 4.4. Let P = (Pi,... ,Pfc) PM? he a collection of d-dimensional lattice poly¬ 
topes. Under which conditions are all (or certain) entries of h*{P) non-negative? Is it 
true that h*(P) > (—1)^'''*(J? What can he said if the polytopes are allowed to he of 
arbitrary dimension? 


In Corollary 4.6 we will show that asymptotically, i.e., if one considers high enough dila¬ 
tions of d-dimensional polytopes, the mixed h*-vector h*(P) always becomes non-negative. 
This suggests that it might be enough to require the lattice polytopes to contain “suffi¬ 
ciently many” interior points. 

Before we provide our main result of this section, we recall the definition of the d**^ Euler¬ 
ian polynomial A^i^z) = X^^=i A(d, k)z^. One, out of several, combinatorial approaches 
to define A,i{z) is the following: 


n>0 


Adjz) 

(1 — z)'^+^' 
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It is known that the Eulerian polynomials Ad{z) have only simple and real roots and all 
roots are negative. The following result, which is a generalization of Theorem 5.1 in [10], 
provides a relation between Eulerian polynomials and mixed ^^-polynomials of lattice 
polytopes. 


Theorem 4.5. Let P 

Then, as r ^ oo, 


= [Pi,..., Pk) be a collection of d-dimensional lattice polytopes. 

JC[k] 


Note that for k = d, (5) implies the result for r = 1. 


Proof. For a d-dimensional lattice polytope P C the Ehrhart polynomial in the usual 
basis takes the form 

Ep{n) = erf(P)n'^-f H-h eo(P) 

where ed{P) = vol(P). For the Ehrhart series we compute 


Y.Ep(n)z’‘ = 

n>0 i=0 n>0 


= 




(1 — zY+^ 

By comparing Erp{n) with Ep{rn), we see that ei{rP) = Pei{P) and hence 


d-l 


Kp{z) = /vol(P)Al,(z) + 5 ^Pei(P)(l-z)''-M,(z). 


i=0 


Now from (10) together with the fact that rP 0 = {0} for all r > 0, we infer 

tkLl = j:±:4 5;(-i)‘-Wvo1(f,)a,(.), 

JC[k] JC[k] 


which shows the claim. 


□ 


For k = 1 this also yields the result of Diaconis and Fulman [10, Theorem 5.1] on the 
asymptotic behavior of the usual h*-polynomial of a lattice polytope. Similar results have 
also been achieved by Brenti and Welker [8] and Beck and Stapledon [4]. 

Before we provide some almost immediate consequences of the previous result, we recall 
that a sequence (oq, Oi,..., a^) of real numbers is called log-concave if > aj_iaj+i 
for all 1 < i < d — 1. The sequence (oq, ai,..., a^) is called unimodal if there exists a 
0 < £ < d such that Oq < ■ ■ ■ < ai > ■ ■ ■ > a^,. 

Corollary 4.6. Let P = (Pi, P 2 ,..., Pk) be a collection of d-dimensional lattice polytopes 
in W^. Then there exists a positive integer R (depending onP) such that for r > R: 

(i) the mixed h*-polynomial h(p{z) has only real roots /d^^^(r) < (3^‘^l{r) < ■ ■ ■ < < 

jd^Yr) < 0 with hmj .^00 /3^'^\r) = pW l < i < d. Here, < • • • < pA'> = 

0 denote the roots of Ad{z). 

(ii) h*{rP) > 0 for 1 < i < d. 

(hi) h*{rP) is log-concave and unimodal. 

Proof, (i) First observe that, by Theorem 4.5, the roots of h*p{z) converge to the roots 
of Ad{z). Moreover, since the roots of Ad{z) are known to be all distinct and negative, 
there exists a positive integer R such that for r > R, h*p{z) has only simple and real 
(negative) roots. Otherwise, since complex roots come up in pairs (if r is a complex root. 
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also its complex conjugate r is a root), the Eulerian polynomial Aii{z) would be forced to 
have a root of multiplicity 2, which yields a contradiction. 


(ii) Using (2) and the fact that ed{P) = vol(P) for a d-dimensional lattice polytope P C 
we recognize the right-hand side in Theorem 4.5 as med(P)Ad(^). By Corollary 2.5, 
med(P) is a non-negative linear combination of mixed volumes of Pi,... ,Pk and thus 
positive. Since the coefficients of the Eulerian polynomials (besides the constant term 


which is 0) are all positive as well. Theorem 4.5 then implies that the sequence 




r>l 


converges to a polynomial, whose coefficients, except for the constant term (which equals 
0), are positive. Hence, there has to exist a positive integer R (depending on P) such 
that for r > P all but the constant coefficient of h*-p{z) are positive. This shows (ii). 

(iii) Using the hrst two parts, we know that there exists a positive integer R such that for 
r > R the polynomial h*p{z) has, except for the constant term, positive coefficients and is 
real-rooted. Theorem 1.2.1 of [7] implies that h*{rP) is log-concave. Since by the choice 
of r, this sequence does not have internal zeros, it is unimodal by [7, Section 2.5]. □ 


Acknowledgment. We thank the referees for very careful reading and their suggestions 
which helped to improve the presentation. 
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